Influence of heavy modes on perturbations in multiple field inflation 
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We investigate linear cosmological perturbations in multiple field inflationary models where some of the 
directions are light while others are heavy (with respect to the Hubble parameter). By integrating out the massive 
degrees of freedom, we determine the multi-dimensional effective theory for the light degrees of freedom and 
give explicitly the propagation matrix that replaces the effective sound speed of the one-dimensional case. We 
^— ^ . then examine in detail the consequences of a sudden turn along the inflationary trajectory, in particular the 

' possible breakdown of the low energy effective theory in case the heavy modes are excited. Resorting to a new 

basis in field space, instead of the usual adiabatic/entropic basis, we study the evolution of the perturbations 
during the turn. In particular, we compute the power spectrum and compare with the result obtained from the 
low energy effective theory. 

oo : 

r 3 1 ' I. INTRODUCTION 

> ■ 

Qj, Although the inflation paradigm has repeatedly confirmed its compatibility with accumulating cosmological data since its 
D ■ emergence more than thirty years ago, the nature of the inflaton field(s) remains elusive. In high energy physics models, it is 
i-^ 1 1 generic to obtain a large number of scalar fields and inflation can be realized if, in this large field space, there exist one or 
several directions along which the effective mass is small with respect to the Hubble parameter H. In addition to these light 
' directions, one would thus find other directions where the effective mass is of the order of the Hubble parameter, or much larger. 
£> . Until recently, the standard approach was just to ignore these massive directions, since the perturbations generated along these 
directions are suppressed. 

Recently, it was realized that the heavy modes can influence the evolution of the light modes if the background trajectory in 
field space is curved. In particular, it was shown in the context of two-field models, that the coupling between the light and 
heavy degrees of freedom can be described as a single-field system with an effective speed of sound that differs from the light 
*n ■ speed when the inflationary trajectory undergoes a turn [Q]], even if the kinetic terms in the original two-field Lagrangian are 
| perfectly canonical. This property was further investigated in J^Ht], as well as in H-p)] . Moreover, when the turn of the trajectory 
, is too sharp, the single field effective description breaks down, as studied in flQIl and full , because the massive degrees of 
freedom become excited during the turn (see also H^-fiH f° r recent works investigating the effects of heavy modes from another 
perspective). 

The goal of the present work is to extend the results of these previous works in several directions. First of all, whereas previous 
studies usually work in the framework of the traditional adiabatic-entropic basis lfl5[|l6ll . we advocate the use of a new basis for 
03 a better intuitive understanding of the behaviour of the background trajectory as well as of the perturbations. This is motivated 
by the consideration that the background velocity, i.e. the adiabatic direction, does not always coincide with the light direction, 
as emphasized in lioll . We thus prefer to use a basis where the distinction between light and heavy modes is manifest. This basis, 
which we call the mass basis, corresponds to the eigenvectors of the effective mass matrix which appears in the Lagrangian 
governing the linear perturbations of the system. In many cases, this basis almost coincides with what we call the potential basis, 
which are the eigenvectors of the matrix whose components are the second derivatives of the potential with respect to the scalar 
fields, i.e. the Hessian of the potential. 

Working in this mass basis, we obtain the effective action governing the dynamics of an arbitrary number of light (i.e. m <C H) 
degrees of freedom in presence of an arbitrary number of heavy (i.e. m ^> H) degrees of freedom. In this generalized framework, 
the effective speed of sound that was obtained previously in the case of a single light degree of freedom, becomes a matrix, whose 
eigenvalues correspond to different effective speeds of sound. 
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This multi-dimensional effective description is valid as long as the inflationary trajectory remains relatively smooth. In some 
models, one could imagine that the potential landscape is sufficiently hilly that the inflationary trajectory has to undergo some 
sharp turns at some particular points. During these sharp turns, one must take into account explicitly the effects of the heavy 
degrees of freedom, as pointed out in ifioll . In the present work, we explore in detail what happens during such a turn, considering 
a two-dimensional description in the two-dimensional plane in field space spanned by the inflationary velocities just before 
and after the turn. We also assume that the trajectory is straight before and after the turn 1 . In order to develop an intuitive 
understanding of the consequences of the turn on the perturbations, we build a simplified analytical model of the turn, which 
contains a parameter controlling the "sharpness" of the turn. For the background, we solve numerically the equations of motion 
in a concrete model, which we compare with our analytical description. 

When the turn is smooth enough, the single-mode effective description is sufficient to compute the power spectrum. By 
contrast, for a sharp turn, it is necessary to go beyond the effective single field description and we obtain an analytical expression 
for the contribution due to massive excitations. Both light and heavy contributions contain oscillations periodic in k. 

The outline of the paper is the following. In the next section, Section [TTJ we introduce our multi-field model and derive the 
equations of motion for the background and for the linear perturbations. Section [Til] is devoted to the derivation of the multi- 
dimensional effective theory when several light directions coexist. In section IIV1 we study the inflationary trajectory during a 
turn. Section [V] deals with the evolution of the perturbations during this turn. We conclude in the final section and give some 
details about our calculations and approximations in the Appendix. 



II. EQUATIONS OF MOTION FOR THE BACKGROUND AND THE PERTURBATIONS 

For simplicity, we restrict ourselves to multi-field models with canonical kinetic terms, although the extension to models with 
non-canonical kinetic terms (following e.g. [20]) should be straightforward. We thus consider N scalar fields (ft 1 , governed by 
the action 

S = j ^x^^-Uud^d^' 1 -V^ 1 )^ , I=1,2,---,N, (1) 

where g is the determinant of the spacetime metric g^ and V(cf) 1 ) is the potential of the scalar fields. We use Einstein's implicit 
summation rule for the scalar field indices. 



A. Background equations of motion 

In a spatially homogeneous and isotropic spacetime, endowed with the metric 

ds 2 = g^dx^dx" = -dt 2 + a 2 (t)8 i jdx i dx :l , (2) 
the evolution of the scale factor a(t) is governed by the Friedmann equations 

H 2 = \ Q<5/ jtftf + V), H= —SutfV = -H 2 e, (3) 

where H = a/ a is the Hubble parameter and a dot denotes a derivative with respect to the cosmic time t. We work in units such 
that Mp = (87rG) -1 ' 2 = 1. The equations of motion for the homogeneous scalar fields are 

ft + 3H<j> 1 + 5 IJ Vj = 0, (4) 

where Vj = dV/d<ft J . 

As we will see later, it is often convenient to make a change of orthonormal basis in field space, so that 

J = <TeL m = l,---,N, (5) 

where the new basis vectors e m satisfy the orthonormality condition 

Sije^ei = 5 mn , (6) 



This differs from a few works I5l^ [l7tjl9ll . where the trajectory follows a circle with constant angular velocity (which we refer to as "constant turn"). 
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as well as <5 mn e^ = 6 IJ . 

In terms of the new components <fi m of the inflationary trajectory, the equations of motion (|4} become 

A<L + 3ff<L + V m = 0, (7) 

where 

D t <f> m = 4>m + Zmn<t>n , Z rnn ■= ^IJ^ m ^ n ■ (8) 

The coefficients Z mn , which are antisymmetric (Z mn = —Z nm ) as a consequence of the orthonormality condition (|6), char- 
acterize the change rate of the new basis with respect to the former field space basis. For example, in a two-dimensional field 
space, the new basis is fully characterized by its angle 9 with respect to the initial basis, and the coefficients Z mn are then directly 
related to 9. 



B. Linear perturbations 

Let us now turn to the linear perturbations of the N scalar fields coupled to gravity. For multi-field models, it is convenient to 
work in the spatially-flat gauge, where the scalar-type perturbative degrees of freedom are encoded in the scalar field perturbations 
Q 1 = 5(f) 1 . The quadratic action for the Q 1 can be obtained by expanding the action 5, given in (fl]i, to quadratic order in the 
perturbations. Using the conformal time r} = J dt/a(t) instead of the cosmic time t, the action for the linear perturbations, 
expressed in terms of the canonically-normalized variables v 1 = aQ 1 , is given by (see e.g. lf2lll ) 

(9) 



S (2) = \\ d V d 3 x [v' T v' - d lV T d % v - a 2 v T (M - H 2 (2 - e)) v] 



where we use a matrix notation, v being the column vector with components v 1 and v T is the corresponding transposed matrix. 
A prime denotes a derivative with respect to the conformal time rj. The matrix Af corresponds to the (squared) mass matrix and 
is given explicitly by 

Mu = Vu + (3 - e) + 1 (Vj4>j + . (10) 



This result has also been extended to models with non-canonical kinetic terms (see e.g. 12011 ). 

If we go to another orthonormal basis e^, as shown in (0, the linear perturbations will be described by new components u m , 
defined as 

?/ = 4X\ m=l,---,N. (11) 

The equation of motion for the u rn can be derived from the quadratic action 

S=\J dvd 3 x u /T u' + u T d 2 u - a 2 u T (m - H 2 (2 - e)) u + 2au' T Zu , (12) 

with 

M:=M + Z 2 , (13) 
where M now denotes the matrix of coefficients M mn (we use the same notation M for simplicity) defined by 

M mn := MijeLd, (14) 



and where Z is the antisymmetric matrix of components Z mn introduced in (|8). In matrix notation, these equations of motion 
are given by 

u" - d 2 u + a 2 (M - H 2 {2 - e) + Z + HZ^j u + 2aZu' = 0. (15) 

So far, our new basis remains arbitrary. The usual approach i s to introduce a basis that consists of unit vectors defined along the 
so-called adiabatic and entropic directions, as first introduced in flal and further developed in lfl6i I22I l23ll . The idea is to use the 
prefered direction in field space defined by the instantaneous background velocity (f) 1 in order to project perturbations parallely 
and perpendicularly to the background velocity (j) 1 , which yields the adiabatic and entropic modes respectively. An advantage of 
this choice is that the comoving curvature perturbation, which is directly relevant for observations, is related with the adiabatic 
mode. Below, we will use of a different basis, which is ideally suited to study perturbations that can be decomposed into light 
modes and heavy modes, as discussed in the next section. 
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III. EFFECTIVE THEORY FOR THE LIGHT MODES 



In the present work, we assume that the scalar field perturbations can be divided into two categories: light modes and heavy 
modes. More precisely, this means that the eigenvalues of the mass matrix M can be organized as follows: 

{to 2 , ■■■ , m\, m 2 L+1 , ■■■ , m 2 N }, (16) 

with 

|m 2 |<ff 2 , (i=l, •••£), to 2 >H 2 , (j = L + 1,---N). (17) 

The corresponding eigenvectors of the mass matrix thus define respectively the light and heavy modes. Note that we do not 
consider here the possibility that modes with a mass of the order of H are present among the perturbations, even if this situation 
can lead to interesting effects as discussed in Il7l [l^l . We also discard the special case where highly unstable modes (with 
to 2 < 0) are present. 



A. Effective action for the light modes 

It is convenient to choose an orthonormal basis that consists of the (normalized) eigenvectors of the mass matrix M, which 
we name the mass basis. In this new basis, the mass matrix is by construction diagonal, i.e. 

M mn = diag{- • ■ , m 2 ;)i , ■ • ■ , m 2 h)j , ■ • • } , (18) 

where the subscript (I) denotes the light modes and (h) the heavy modes. The corresponding perturbation modes u m in (fTTT > can 
also be divided into light modes u^ m and heavy modes U(h) m , i- e - 



u (i) 
u {h) 



(19) 



As a consequence, the antisymmetric matrix Z can be decomposed as 

Z d) Z M 



Z = 71 • (20) 



~ Z (c) Z (h) 

In order to determine the effective action describing the dynamics of the light modes, we follow the strategy of JH 0]. By 
neglecting the kinetic terms of the heavy modes in (fT2] i. the variation of the action with respect to Ufh\ yields the constraints 

- Au (h} + \ z f c) u' {l) + O t u {1) = 0, (21) 

with 

A := M (h) - Zf c) Z {c) -H 2 (2-e)-^, (22) 

O := HZ( C ) + Z( c ) — Z(/)Z( C ) - Z( c )Z(/i). (23) 

Plugging the solutions of the constraints, 

u (h) = A- 1 (^(Vffl + ° Tu w) > W 

back into the action (fT2l gives, after some straightforward manipulations, the effective action for the light modes which can be 
expressed in the form 



Srff[«(o] = ^ / dr]d 3 x u'^Ku' ([) + uj^uy) - a 2 u^ (M - H 2 (2 - e)) uy) + 2au'^ ) Zu(i) 



(25) 



The coefficients of the above action are explicitly given by the expressions 

7C = l + 4Si, (26) 

M = M (l) + Z 2 } - Z (C) Z ( T C) - H 2 (Hi + H 2 + + S 3 ) + ^ [aH (2Si + H 2 + S^)]' , (27) 

Z = Z(;) + H (H 2 - H^) , (28) 
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with 

Hi := Z( C )A Z( c ), H 2 := — Z (c )A ^Z( c ) — Z(;)Z( C ) — Z^Z^J , (29) 

H 3 : = ^2 (%) - z (i) z (c) ~ Z(c)Z(h)) A -1 (Z (c) - Z (i) Z (c) - Z (c) Z (/l) ) , (30) 

Note that all Sj's are suppressed by A -1 , which is defined in d22t . In the regime where we expect the effective theory applies, 
the eigenvalues of JVf^) should be much larger than the elements of the other matrices that appear in A, and one can make 
an expansion in the small parameter , where m\ corresponds to the typical value of the eigenvalues of Mqa. Thus one 
can write A -1 = M7j}(l + 0(m^ )) and, at leading order, replace A^ 1 by MT^; in Hj's, which make the expressions much 
simpler. For later convenience, we introduce the notation 

A = Z (c)M^Z( c) , 

and thus write 

St = A + 0(m~ 4 ). 

The action (l25l l is one of our main results in this work. The effects of the heavy modes on the light modes are effectively 
encoded in Z c and the dimensionless matrices Hi, E 2 and H3. One notices that the influence of the massive modes multi-field 
manifests itself only when the background trajectory is turning in field space, i.e. when the mass basis undergoes a rotation, so 
that the coefficients Z nm are nonzero. 

B. Discussion on the sound speeds 

From ( I25I 1. the light modes are propagating with effective speeds given by the eigenvalues of 

/C- 1 = (l+4Hi)- 1 =l-4A + G(m- 4 ). (31) 
In the following we discuss some particular examples. 



1. Single light mode 

When one assumes that there is only one light direction, the above matrices become numbers. One then finds that the light 
mode is characterized by a propagation speed given, up to order in m^ 2 , by 

N Z 2 

^^ 4 E^ + K- 4 ) (32) 

j=2 ( h )i 

where the sum is earned out over all massive modes. 

In the particular case where there is only one heavy mode as well, one immediately recovers the results of H. Indeed, since 
Z( c ) coincides with —6, where 9 denotes the angle of the adiabatic direction in field space which is commonly used in the studies 
of two-field models, one finds that the matrix JC reduces to a number given by 

1 _ 40 2 

~% ~ 1 + ml - 02 _ (2 - e )H* + g ' (33) 
where k is the wavenumber associated with the Fourier modes of the perturbations. 



2. Several light modes 

When several light modes are involved, one can obtain an anisotropic propagation of the perturbations, similarly to what can 
be found with generalized multi-field Lagrangians with non-standard kinetic terms (see 1I20I0 . The corresponding sound speeds 
can be obtained by diagonalizing the matrix A, whose explicit elements now read 

A«= £ ^ (34) 

mi 

h=L+l h 
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where the index h is summed over the heavy field indices. 

Denoting Aj the eigenvalues of the matrix A, one finds that the associated sound speeds are simply given by 

4 = l-4A l + 0(m,; 4 ), (35) 

which is the direct generalization of (|32V 

Note that if J2i ul Zih = for some non-vanishing (light) vector, then Ay becomes degenerate, with u l as eigenvector 
associated with the eigenvalue A = 0. This means that the propagation speed along the u l direction is the speed of light as usual. 
This is the case if a given light mode is not coupled to any of the heavy modes (i.e. if Zui = for all values of h). Note also 
that, if there are as many light modes and heavy modes, Zi c \ is a square matrix and dct(A) = dct(Z^) 2 / det(Af( ft )). If 
is degenerate, then A = is an eigenvalue of A. 



IV. SUDDEN TURN ALONG THE INFLATIONARY TRAJECTORY 



In the present section, we focus on the possibility that the inflationary trajectory undergoes a turn due to the sudden change of 
orientation of the massive directions, and correspondingly of the light directions. It is a reasonable approximation to restrict the 
analysis to a two-field description in the two-dimensional plane spanned by the velocities before the turn and after the turn. We 
will also assume that there is enough time between two sudden changes of direction that the inflationary velocity can relax to its 
natural direction between two turns. 

During a turn, the trajectory deviates from the bottom of the potential valley because of the centrifugal force, as illustrated in 
FiglH If the turn is sufficiently soft then the trajectory smoothly relaxes to the minimum of the potential after the turn. However, 
if the turn is sharp, the trajectory strongly deviates from the bottom of the valley and then relaxes via oscillations around the 
"minimal" trajectory. In the latter case, the heavy fields become excited and the amplitude of the oscillations depends on the 
"sharpness" (defined precisely below) of the turn. 



AO 



(T 



FIG. 1 : Schematic representation of the background trajectory in field space dur- 
ing a single turn. The dashed line represents the light direction of the potential 
(i.e. the bottom of the valley). If the turn is "soft", the trajectory deviates slightly 
from the bottom and remains smooth (blue line). If the turn is sharp, the trajectory 
deviates from the light valley significantly and oscillates after the turn (red line). 

When the relaxation is oscillatory, the direction of the inflationary velocity strongly varies during the process. This indicates 
that choosing the adiabatic-entropic basis to describe the inflationary evolution might not be the most appropriate as this basis 
undergoes wild oscillations during the turning phase. This is why we will use the potential basis, which does not oscillate during 
this process, to study in detail the background trajectory during and just after the turn. The evolution of the linear perturbations 
will be investigated in the next section. 



A. Kinematic basis 



Before using the potential basis in the next subsection, let us first describe the background evolution in terms of the usual 
adiabatic-entropic decomposition. The adiabatic direction is represented by the unit vector 



n =y, &= ySutfp, (36) 



while the entropic direction is along the unit vector s 1 , which is orthogonal to n 1 . We will call the basis {n 1 , s 1 } the "kinematic 
basis", since n 1 is pointing in the velocity direction. 

Starting from the background equation of motion (0) in the original basis (or from (|7) in an arbitrary orthonormal basis), the 
projection onto the adiabatic (velocity) direction, whose components are denoted n 1 (or n m ), yields 

a + 3Ha + V. a = 0, V >a = rJVj . (37) 

Since n 1 is a unit vector, its time variation is orthogonal to n 1 , and therefore proportional to s , so that one can write 

h 1 := 6s 1 (38) 

where the coefficient of proportionality, 9, corresponds to the time derivative of the angle 9 of n 1 with respect to the initial field 
basis. 

One can also derive a second order equation of motion for 9, by using the projection of the background equation (|4} along the 
entropic direction, which reads 

&6 + V.s =0, V s = s'Vj . (39) 

The time derivative of this equation yields 

&6 + a9 + &V as -9V a =0, V as = n 1 s J V,u , (40) 
where we have used s 1 = —On . Dividing by a and using the adiabatic equation (|37| | to eliminate V a , we finally obtain 

e+(sH + 2^je + v i(TS = o. (4i) 

B. Potential basis 

We now restrict ourselves to two-field models, working in an arbitrary orthonormal basis where the components of the adiabatic 
and entropic unit vectors are denoted n m and s m , respectively. In this case, the entropic subspace is one-dimensional, along s m 
which can be expressed as 

s m = -e mn n n , (42) 
where e mn is the fully antisymmetric two-dimensional tensor (with en = 1). This implies in particular 

V as = n m s n V mn = n x n 2 {V 22 - V n ) + (n? - A) V, 12 , (43) 
where the components of the adiabatic unit vector can be written in terms of an angle ip 

{m, n 2 } = {cos ip, sin tp} . (44) 

So far, we have been working with an arbitrary orthonormal basis. It is now convenient, when there exists a mass hierarchy 
between the various directions in field space as assumed in the present work, to choose explicitly the "potential" basis, i.e. the 
eigenvectors of the Hessian matrix of the potential 2 . In the potential basis, we thus have 

V mn = diag{mf ,m|}, with m; < H < m h . (45) 

The angle ip introduced in d44l i corresponds to the angle of the background velocity with respect to the "potential basis". The 
angle 9 is thus the sum 

9 = 4> + 8 pi (46) 
where 6 P is the angle between the potential basis and the original field basis. 



2 Strictly speaking, the potential basis is different from the mass basis. However, when there is a large mass hierarchy, the potential basis and the mass basis 
almost coincide in general because, in the heavy subspace, the first term on the right hand side of dlOl is much bigger than the other terms. 



Substituting this decomposition into d4TT >. one obtains 



ip + 3H 1 



2a 
3Ha 



ip+- (ml - mf) sin(2i/0 



3H 1 



(47) 



i.e. an equation of motion for ip, where the source term depends on the angle 6 p . The interest of this equation is that, in general, 
one expects the behaviour of 6 p to be rather smooth whereas the angles ip and 6 will oscillate wildly if the inflationary trajectory 
is strongly displaced from the local minimum of the potential. With the above equation, one can simply interpret the evolution 
of tp as a response to the change in the potential landscape. 

In order to obtain analytical solutions for (147) , we will assume that the inflationary trajectory satisfies the slow-roll approxi- 
mation, so that a / (Ha) can be neglected. We also assume that the angle ip is sufficiently small so that the sinus can be replaced 
by its argument. With these approximations, (l47t reduces to 



ip + ZHiJ; 



m\ip = -9 p - 3H9 P 



(48) 



where one recognizes the familiar equation of motion of a damped harmonic oscillator with a source term that depends on 9 p (t). 
Assuming that H and remain approximately constant during the turn, (|48] can be formally solved as 



dt'G(t,t') 6 p (t') + 3H6 p (t') 



where the retarded Green's function is given by 

G(t,t') = Q(t-t') ShLiu;{t - t ' )) e-i H ^'\ 

UJ 

where 6 is the Heaviside distribution. 



(49) 



(50) 



C. Modelling the turn 

If the potential is sufficiently regular, the evolution of the angle 6 p (t) during the turn can be described by the expression 

e p (t) = MS(t-to), (51) 

where AO is the global variation of the angle during the turn, to represents the turning time, and S(t) is some function which 
interpolates smoothly from to 1 . 

The precise shape of the function S(t) depends on the details of the potential, but one can try to understand the generic features 
of the behaviour of %p by introducing a simple Gaussian shape for the time derivative of S(t), so that 

9 p (t) = A6-^e-^ 2t ' 2 . (52) 
V 2tt 



We have fixed to = for simplicity and the parameter /j, characterizes the duration of the turn: A£ tm -n ~ A* -1 - We thus have 
three characteristic energy scales: H, w(i=a irih 3> H) and /i. 
Plugging ( T52b into d49l , we get 



AO 9H 2 

m = — 



^l + ^ e -*« 5? [e^ 2 / 2 e-^erfc (-^) 



(53) 



where 5i denotes the real part of the argument and erfc(z) = 1 — crf(z) is the complementary error function. We have also 
introduced the parameters 



Note that, since lo ^> H, we have a w 3H/ (2lu) <C 1 and ip w iui/fj,. 

It is now useful to distinguish two distinct qualitative behaviours depending on the value of the ratio ju/w: 
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FIG. 2: Evolution of the angle ip during a sharp turn as a function of the number of e-folds, 
for various values of fi/io (and rrih /H = 10). The colored continuous lines correspond to the 
numerical integration of Eq.(l47t with constant H , while the black dashed lines are obtained 
from the analytical approximation Eq, d55k 



1. Sharp turn (n/ui > 1) 

In this case, ip(t) can be approximated by 



AO _j£ ( at\ s Hf ( 3Hlj\ 
A= e~* m cos cot - a 



1>V) « e ^crfc ^-^=J ^ cos [cot- a- j . (55) 

The evolution of ip(t) is illustrated in Figj2] Long before the turn, — fit ^s> 1, erfc(— fit/y/2) — > and thus ijj — > 0. Soon after the 
turn, erfc(— /Mt/v2) — > 2, and ^ starts to oscillate with the damping factor e~? Ht . Note that the initial amplitude of oscillation 

is approximately ~ e *i>? A8 < A8. 

In the limit \i ^> cj, the expression dB31 l further simplifies into 



tp(t) = -Aee-^* cos (cut- a) Q(t) , (56) 

which can also be obtained directly from (|49l by assuming an instantaneous source term, i.e. S(t) = 0(f) and therefore 
9 p {t) = AO Sit). For an instantaneous turn, the amplitude of ijj(t) is maximal at the turn and then oscillates with an amplitude 
damped by Hubble friction. 



2. Soft turn ( ft / 'id <C 1) 
In this limit, before and during the turn, (IB3T l can be approximated by 

m^^e-^Lt-Z*), (57) 

V27T UJ Z \ fX J 

which shows that there is no oscillating response during a soft turn, as illustrated in Figj3] The largest value of \ip\ occurs around 
the time 



3H / 4u 2 

'*-2^' 7:= V 1+ 9^- 1 ' 



which is before the turning point to = and the corresponding value of ip is 



3-ff/z 



Vmin = « ~Ag - (2 + 7 ) e"^ 7 . (59) 
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FIG. 3: Evolution of ip during a soft turn (for uiyJH = 10). The colored lines correspond to 
the numerical integration of Eq. d47t with constant H, while the black dashed lines are given 
by the analytical approximation EqJ57t. 



Since uj S> fi, H, the above amplitude is much smaller than AO. 

Long after the turn, when fit > u>/ fi S> possesses an oscillating component but with a tiny (and thus negligible) amplitude 
since w/n> 1: 



AO 



e 2 ** e 



cos [cot — a 



3Hlu 



1 



2ir 



2 1 
fit 



To summarize, in a soft turn, the background trajectory is almost the same as the light direction of the mass basis, i.e. 
kinematic basis is tightly following the mass basis. There is no explicit oscillation of the trajectory. Moreover, the "softer" 
turn is, the closer the background trajectory is to the direction of mass basis. 



(60) 

the 
the 



D. Numerical study 

In order to confirm the qualitative conclusions of our analytical approximations, we have studied numerically a concrete two- 
field inflationary model with a potential that contains a valley with a turn. Our potential for the scalar fields (f> and \ 1S given 
by 

V(<t>,x) = \m 2 cos 2 (^pj [ X -(0-0o)tanS] 2 + i™^ 2 , (6 1) 

where A8 is a constant parameter corresponding to the net angle variation due to the turn. The function S is defined by 

AO 



TV 



arctan[s(0 — 0o)] , (62) 



where s is a constant parameter that controls the sharpness of the turn. By choosing M 3> m$, as illustrated in Fig. @] this 
potential yields a steep-sided valley, characterized by 

X = (4>- <MtanS, (63) 

which undergoes a turn around 4> = 0n, X = 0. 

To make easier the comparison with our analytical model of the turn, it is useful to obtain a rough relation between the 
parameters s and ft. This can be done by approximating the function arctan with the error function erf, more precisely 
arctan(a;) rj 7rerf (x)/2, where the proportionality coefficient is chosen so that the asymptotic limits coincide. This yields 

AO AO f 3 ^-^ 2 

S « —evi(s(cf> - cf> )) = -= / e-y-dy. (64) 

2 VTT./n 
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FIG. 4: Contour plot of the valley corresponding to the potential 
Eq. d61t . with the numerical parameters given in Eq. d68t . 



Noting that when the inflationary trajectory evolves along the valley, the potential basis angle 6 P coincides with 5, one obtains 
the following approximation for ( 



3« A0- 



-s 2 (<j>-<j> ) 2 



(65) 



Comparing this expression with Eq. (152l i. we find the relation 

jj, w splits . 



(66) 



if <f> is expressed as a linear function of t (note that </> > 0) Furthermore, if the last term in Eq. (l6ll dominates the potential and 
the slow-roll condition is satisfied for <j>, /i can be expressed as 



2y/3 



-m^s . 



(67) 



We have integrated numerially the background equations of motion for the scalar fields and x m tne potential (|6*TT i with the 
following parameters 



M = 1.0xl(n 4 , m = l.OxlCT 7 , 0o = -lOth/rj, A6» = 7r/10. 



(68) 



expressed in Planck mass units (i.e. Mpi = 1) for various values of the sharpness parameter s, in particular for s = 100\/3, 
s = 500\/3 and s = 1000V3, which correspond to the values fj,/u) = 0.2, (j,/oj = 1 and (j,/u) = 2, respectively, according to 
Eq. (ISTl i. In all cases, we start our integration at the position (</>, \) = ( — 2 — lOOv^, 2 tan(?r/20)) with vanishing velocity (the 
initial velocity is in fact irrelevant since inflation quickly reaches the slow-roll regime, long before the turn). 

The background trajectories in field space around the turning point for three different values of s, are shown in Fig. [5] The 
turn becomes sharper for higher values of s. For s = 100V3, the trajectory is smooth, but oscillations are clearly visible for 
s = 500^ and s = 1000\/3. This is even more conspicuous when one considers the inflationary velocity a as shown in Fig. [6] 
For the three inflationary trajectories, we have checked numerically that the slow-roll parameter e remains very small (even if 
some oscillations are observed for s = 500^ and s = 1000>/3). 

Along the trajectory obtained numerically one can compute the evolution of the eigenvalues of the mass matrix which can be 
denoted mf and as usual. One can check that the highest eigenvalue nth is of the order of M, as expected. For all the three 
trajectories considered, the hierarchy mf < H 2 < m\ remains valid during the turn. However, because of the specific shape 
of the potential around the turn, the simplifying assumption that m/j is constant is not strictly valid: we observe the reduction 
(of about 2.5 %) of rrih during the turn. Another subtlety is that mf becomes negative just at the moment of the turn, with an 
amplitude that increases with the sharpness of the turn. 

Let us finally discuss the evolution of the kinematic basis, chracterized by the angle 9, with respect to that of the mass basis, 
which almost coincides with the potential basis 3 . For a smooth turn (s = IOOa/3), the velocity has the time to adjust to the 



1 We have checked explicitly that the potential and mass bases almost coincide: during and after the turn, we find that 



- Sm)/ is at most of order 10 
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FIG. 5: Background trajectories around the turn for different 
values of the sharpness parameter s. The field (f> increases 
with time. 




AN 



FIG. 6: Evolution of a around the turn for different values of 
s. 



rotation of the mass basis and the two angles evolves together as illustrated in Fig. [7] By contrast, for a sharp turn, the velocity 
follows the light direction with some initial delay and then oscillates, as illustrated in Fig. [8] for IOOOa/3. 

For comparison with our previous discussions, we show the evolution of ip = 6 — m in Figs. [9] and [10] for the cases with 
a sharp turn and a soft turn, respectively. One sees clearly that the evolution of ip is very similar to that of our approximate 
analytical solutions. 
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FIG. 7: Evolution of 6 and m around the turn for s 
100^3. 



FIG. 8: Evolution of 
1000^. 



and 9 m around the turn for s 




FIG. 9: Evolution of ip during a sharp turn for different val- 
ues of s, corresponding to h/uj — 1, 2 and 5, respectively, 
according to Eq. doTt . 



FIG. 10: Evolution of ip during a soft turn for different values 
of s, corresponding to [i/uj = 0.1, 0.15 and 0.2, respectively. 
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V. EVOLUTION OF THE PERTURBATIONS DURING A TURN 

After the study of the background evolution, we now investigate the evolution of perturbations during the turn in order to see 
whether specific signatures will be imprinted in the final power spectrum. 

A. Two-dimensional analysis 

I. Adiabatic-entropic basis versus mass basis 

Let us start with the familiar approach where perturbations are described in the adiabatic-entropic decomposition, i.e. the 
kinematic basis. The equations of motion for the adiabatic and entropic modes read lfl5ll . in Fourier space, 

2" (zO'Y 

u" + k 2 u a u a = 2e'u'+2± '-u s , (69) 

z z 

u" s + k 2 u s + a 2 (v ss - 9 2 - H 2 (2 - e)) u s = -26'u' a + 2—6'u„, (70) 

where z = a\[2e and 9 is, as previously, the angle of the adiabatic direction in field space. As we saw in Sec II V Cl during a sharp 
turn, the trajectory deviates from the bottom of the valley, and then rapidly oscillate around the curve that follows the bottom of 
the valley. As a consequence, the angle 6 undergoes strong oscillations just after the turn, which makes the analysis of the above 
system of equations somewhat complicated. In this respect, the situation becomes much simpler when working in the mass basis, 
because the evolution of the corresponding 6 m evolves much more smoothly as we have seen earlier. 

It is thus interesting to study the equations of motion for the perturbations described in the mass basis. Using ( fT5t , one 
immediately finds 

+ 2fi>h. (71) 

- 20>;, (72) 

The smoother behaviour of 8 m will enable us to obtain approximate analytical solutions for the system (I7l1i-(l72l. as we will see 
below 4 . 

Note that the descriptions in the kinematic basis and in the mass basis are simply related by a rotation of angle rj) = 9 — 9 m , 
so that, at any time, the adiabatic mode, for instance, is given by 

u a = cos ip ut + sin ip Uh , (73) 

in terms of the light and heavy modes. Moreover, long after the turn, the angle ip relaxes to zero if we assume that there is 
no subsequent turn. Therefore, the light mode and the adiabatic mode will coincide in the asymptotic limit, which allow us to 
compute the prediction for the adiabatic power spectrum either by using the adiabatic mode, or the light mode at late time. 

2. Three ranges of scales 

Before deriving approximate solutions for the system in the next subsection, it is useful to identify the main relevant scales. 
In standard inflation, it is convenient to distinguish, at a given time, sub-Hubble and super-Hubble perturbations. In the present 
context, it will be useful to distinguish three intervals of wavelengths, or wave numbers, because the heavy mode introduces the 
new characteristic mass scale m^, in addition to the usual cosmological scale H. 

Our system (l7lTi-(l72l thus appears to have two "horizons": one is the usual Hubble horizon with comoving radius (aH)^ 1 , 
the other is the "heavy mass horizon" with comoving radius (am^ 1 . Since mu ^> H, the mass horizon is always deep inside 
the Hubble horizon. During their evolution, the modes will therefore experience two "horizon-exit" processes, as illustrated in 
FigJTTJ first the "mass-horizon" exit, and then the usual Hubble exit. Denoting a* the scale factor at the time of the turn, one can 
thus distinguish three categories of modes: 



" + [k z + a z mf 



k 2 + a 2 m 2 - e' 2 



a 
a 

a 



4 Note that 8 m (and also 9 P ) has an oscillatory component, because the trajectory oscillates around the bottom of the valley, but this oscillatory component is in 
general subdominant with respect to the smooth component. 
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FIG. 11: Evolution of the comoving Hubble radius (aH)^ 1 and of the "heavy mass 
horizon" (arrih) ~ 1 • The Fourier modes, characterized by a fixed comoving wavenum- 
ber k, exit the heavy mass horizon first and the Hubble radius later. 



• subhorizon modes fc/a* S> rrih (region I of FigJTTb: these modes are deep inside both the Hubble horizon and the "mass 
horizon", when the turn occurs. 

• intermediate modes rrih 3> fc/a* H (region II of FigJTTb: these modes are inside the Hubble horizon but outside the 
"mass horizon". 

• superhorizon modes H ^> fc/a* (region III of FigJTTb: these modes are already outside the Hubble radius at the time of 
the turn. 

In the next subsection, we investigate in detail how the power spectrum is affected by the turn in the context of our simplified 
analytical model of Section HV CI considering successively the three intervals of wavenumbers discussed above. 



B. Approximate analytical solutions 

The system of equations d7Tl)-(l72b requires in general a numerical treatment. However, it is useful to derive approximate 
analytical expressions in order to gain an intuitive understanding of the qualitative behavior of the system. In order to do so, we 
use a perturbative approach, by noting that the coupling between light and heavy modes is effective for a short duration only, 
since 9' rn is significant only during the short time interval corresponding to the turn. 

In absence of the turn, i.e. for 6' m = 0, the light and heavy modes satisfy the decoupled equations of motion 

< (0 ) + [k 2 + a 2 mf - ^ u m = 0, < (0) + (k 2 + a 2 m 2 h - u h{0) = . (74) 

In the following, for simplicity, we will assume that m; can be neglected, rrih is constant and the background evolution is 
approximately de Sitter. Using the usual Bunch-Davies initial conditions, the relevant solutions for the above decoupled system 
are then 



j — ikr) 



U m(^) = - m [i--) 1 (75) 



and 



u H0) (V, k) = £Le-W^H® {-H , ^V^-I = 77"9 i>>1 ' (?6) 

where i?W denotes a Hankel function of the first kind. 

The effect of the coupling between the light and heavy modes can be taken into account by solving (f7TTi-(f72]i perturbatively 
via Green's function method. This means that a solution of the first equation in d74l with a source term S on the right hand side 
can be written as the following linear combination of positive and negative frequency modes: 

ui(r),k) ~ u m (ri,k) + C(r), k) u m (rj, k) + D(rj, k) < (0 )(r?, k), (77) 
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with 

C(r),k)=i P d V 'u^ 0) (ri',k)S(r]',k), D( V ,k) = -i J dr,'u m (r]',k)S(ri',k). (78) 

The source term must be sufficiently small for the perturbative treatment to apply. In the following, we will use our simplified 
analytical description of the turn ((52), assuming that 0, p and m coincide, and treat AO to be a small parameter. On the right hand 
side of the equation of motion for the light mode, it is useful to distinguish the light contribution to the source term 

Si = e'lu m , (79) 

which is proportional to (AO) 2 , and the heavy contribution, 

S h = O'L u h(o) + 20' m u h{0) , (80) 

which is linear in the small parameter A9. 

We are ultimately interested in the power spectrum of the light mode on super-Hubble scales since the light mode coincides 
with the adiabatic mode after the relaxation phase following the turn. Because the light and heavy modes are initially statistically 
independent, the final power spectrum is the sum of two contributions: 

V(k) = V {l) (k)+r {h) (k) = lim -|L (| u fVfc)| 2 + h ( %,fc)| 2 ) . (81) 

The first contribution is obtained by putting the heavy mode to zero initially while taking Bunch-Davies initial conditions for the 
light mode. The second contribution comes from initial conditions where the light mode is set to zero while the heavy mode is 
in its Bunch-Davies vacuum. Using d77i i, each contribution can be written as 



V {l) (k) = 1 + Cr<0 -Z>W V V {h) (k) = lim J^\ u \ h \r,,k)\ 

/c|t|<§C1 Z7T a 



c (h) _ D (h) 



(82) 



where C^- 1 '^ and D^ l ' h ^ denote the respective limits of the expressions (178) when 77 — > 0, and Vq is the power spectrum arising 
from Ui(o), which is here scale-invariant. 

For the the spectrum due to the heavy modes, the corresponding coefficients and are obtained by plugging Sh into 
(ED. This yields 



V [h) {k) = \I h {k)yV Q , I h (k) = 2ildr)'$l[u m (ri',k)]S h (r)',k), (83) 

which is proportional to (AO) 2 . The light contribution to the spectrum contains a zeroth order term arising from Uhq} and 

perturbative corrections which start at quadratic order in AO, since the coefficients c[ and D± \ obtained by plugging Si into 
(|78l , are already proportional to (AO) 2 . There is no term linear in AO, but there is another term in (AO) 2 , which comes from the 
equation of motion for the heavy mode. Indeed, even if the heavy mode is initially zero, the source term on the right hand side 
of d72l will generate a heavy mode proportional to AO. This heavy mode will in turn produce, via the coupling on the right hand 
side of (l7Tl i. a new contribution to the light mode which is proportional to (AO) 2 . This new term is more complicated to evaluate 
as it arises from a double integration but is of the same order of the other contributions and must therefore be include! d. 

In summary, one finds that the leading order perturbation to the power spectrum, proportional to (AO) 2 is the sum of three 
contributions, given by 

AV 

Ti + T h + Tih (84) 



with 



V 



Fi = 21St[Ii] h(k) = 2ij dr)'R[u m (r,',k)]Si(r,',k), (85) 



J 7 h = \h(k)\ 2 h(k) = 2i f dT/»[u l(0 )(»/,fc)]Sh(»/>*). ( 86 ) 
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and 



ii, 



2R[J Jh ], Ji h {k) 



drnM 0) M C (m) «r fo) + 2 C (m)«r fa) 



(0) 



'(0) 



'/I 



x / d m uf* ( m ) C M uY 1 {m) + 20'm im) uT {m) 



(0) 



,'(0) 



'/i 



x / fe) C ( V2 ) u\ 0) ( m ) + 28' m (772) «; (0) (r? 2 ) 



(87) 



In the following, we give some analytical expressions for the integrals and Ih, expressed in terms of A8, of the dimensionless 
parameters 



A- 



H ' 



(88) 



as well as v, already defined in ( |76| |. These approximate expressions for the integrals // and Ih are obtained by writing 8', which 
appears in the source terms Si and Sh, in the form 

8' m = a8 m = a.A^e-^ 2 "-'-) 2 ^"-'-) . (89) 

V 27T 

The integrals and also contain the uncoupled mode functions Uum and Uh(o)> given in d75l l and d76b respectively. In order 
to perform analytically the integrals // and 1^, it is convenient to consider separately the three ranges of scales, where simpler 
expressions for Uuq\ and uuq) can be used. 



1. Scales k/(a*mh) 1 (orx, 3> v) 

In this case, modes are inside the "massive horizon" (and thus also inside the Hubble horizon), which corresponds to region I 
in figure [TT] Both light and heavy unperturbed mode functions are approximated by 

e — ikr) 

Uj(0) ~ u h(0) ~ —7=- ( 9 °) 
V Ik 



Assuming ^> ft 2 ^> 1, in order to get simple expressions (see the appendix for details), one finds that the dominant contribu- 
tions to the light and heavy integrals are given by 

iz(fc) « i(A8) 2 — (ev? + e -^-^e 1 - il2+tx A , I h (k)^A8e^, (91) 
4y7rx* V 



where one notices that Ih does not depend on v. Since the real part of /; is suppressed by the exponential factor cxp(l — jj 2 ), the 
contribution Ti becomes negligible while 

T h « (AS) 2 . (92) 



2. Scales H <C fc/a* <C m h (or 1 < x, < v) 

In this case, modes are outside the "mass horizon" but still inside the Hubble radius. In this regime, the heavy mode function 
(|76*| | can be replaced by 

-ii/(ln(2i/)-l) 

Uh(0) « /= e- 1 "^^, (93) 

where we have used the approximation ll24ll 

^"r^ijCff > V>>X>Q - (94) 
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as well as 

r(iu + l)^e-i wv+ i hl ^e i d +u(lnu - 1) \ v > 1 (95) 

for large positive v. Using the above expression, we obtained explicit analytical expressions for // and Ih(%*, £t, v), which are 
given in the appendix. As they are rather complicated, we discuss here, in the main text, only the regime x„ <C jl, in which the 
dominant contributions to the integrals reduce to the simple expressions: 

Ii ps -7^7= J 1 ' e lx " (x* + i) (x* cosx* — sinx„) , (96) 



A9y^ 

x 3/2 



e'^ r e ll/ln ^ (x*cosx* -sinx*). (97) 



One thus finds that // is an oscillatory function of x* (note that in the general case discussed in the appendix, I\ also contains a 
non-oscillatory contribution). The integral Ih contains an additional oscillation periodic in v In x* but this will disappear in the 
power spectrum which depends only on the modulus of Ih- Concerning the amplitude of the light and heavy contributions, one 
can see that for large values of x* (if v is sufficiently large), they behave as 

~ (A6>) 2 — , I h ~Me~&tf^. (98) 

X if; y X if; 

This implies that both amplitudes decrease for higher values of k, the suppression being more pronounced for Let us also 
notice that the heavy contribution is strongly suppressed when fi <C v, which corresponds to a soft turn. 
Substituting d96i l and d97l i into (I85H861 I , one finds the following perturbations to the power spectrum: 



Ti ps — (Ad) 2 !2 (x» since, + cos a?*) (x» cos a;* — sin a;*) , Th & (A0) z — ~e v ^ (x* cosa;* — sina;*) z .(99) 



— — (t.~ sin t. -I- msx.^ fi. mi. — sin t..A . 7-Y ( A#"l 2 — p. v ^ 



r 



The deformed power spectrum is illustrated in FigJT2lwith three different choices of jl, corresponding to increasing sharpness. 
The parameter A9 = tt/30 is sufficiently small to find an excellent agreement between the perturbative approach and the full 
numerical result. As shown on the right panels of Fig{l2] the deformation of the power spectrum is dominated by the light 
contribution T\ for a soft turn, but the two other contributions Tk and T\h become significant for sharp turns. Moreover, 
whereas Th is necessarily positive, the terms Ti and Tih can be positive or negative. 



3. Scales k/{a*H) < 1 (or x* < 1) 

For these modes, the turn occurs when they are already far outside the Hubble radius. One can thus obtain analytical ex- 
pressions by resorting to an expansion with respect to the small parameter x*, as discussed in the appendix. This leads to the 
integrals 

h = (Ad) 2 -^= (3 + ie^xl), (100) 



18V^ 

Ih = -Ag {2V + : U) e ^ + H 1+ ^) H9t) . (101) 

The purely light and heavy contributions to the power spectrum perturbations are thus given by 

(A8) 2 (AO) 2 . , N 2=i£i , 

T X = { -—Lfl , T h = h-^ 4^ 2 + 9 x\ . (102) 

The first contribution is fc-independent term, whereas the second contribution is proportional to k 3 , which is suppressed in the 
super-Hubble domain considered here. Note that the latter term, due to the heavy modes, is further suppressed for soft turns via 
the exponential factor cxp(— to 2 //i 2 ). 



C. Comparison with the effective theory 

It is instructive to compare the results that we have just obtained in the full two-dimensional analysis with the predictions of 
the effective field theory presented in Section|III] In the present case, the effective theory contains a single mode and the equation 
of motion for the corresponding canonical variable v = it; /c s is given by 

c 2 k 2 + a 2 m 2 s - — )v = 0, (103) 
a J 
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FIG. 12: Power spectra for the values jl — 12, 60, 120, with the other fixed parameters fixed to the values v = 60 and A6 = 7r/30. 
On the left panels are plotted the fully numerical result (red), the result from the perturbative approach (green) and the result from the 
effective theory approached (dashed blue). On the rights panels are plotted the three contributions Ti (blue), Th (red) and Tin (green) to the 
perturbation of the power spectrum, i.e. the sum of these three curves (plus one) coincides with the green curves of the left panels. 



with the effective speed of sound 



1 + 4- 



to? = to? — at. 



H A (2 - e) 



,,2 



(104) 



and the effective mass 



'elf 



2„4 



rata 



a \ ra t a 



2 2 m r n 



(105) 



In the above equations, 6 m denotes, as usual, the angle of the mass basis with respect to the original basis. 

As before, we assume for simplicity that H and rat are constant, and m; = 0. In the regime where the effective theory is 
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expected to be valid, i.e. m\ 3> 9 2 n , we may expand both c s and m e g with respect to large m/j to obtain 



c 



J i-^i + ^K 4 ). (106) 



'ft 



y'2 



m off * + i i^ 2 H 2 9' 2 + 4^ + 12a^CC - K? - 2 CC) + O (m* 4 ) . (107) 

In the decoupling limit where rrih oo, one finds c 2 . — s- 1 and m^ ff — > — # 2 n , and the equation of motion d 1 03t is then equivalent 
to the equation of motion for the light mode (|7T1 i with the heavy mode Uh put to zero. 

In order to obtain analytical solutions, we can solve (1 1 031 > perturbatively via Green's function method, as we did before. For 
this purpose, one rewrites (11 03t in the form 



v » + (k 2 -?p)v=[(l-c 2 s )k 2 -a 2 m 2 eS ] v, 



(108) 



where the right-hand side is treated as a source term. The unperturbed solution U( ), i.e. the solution of the homogeneous equation 
(with the right-hand side put to zero), is the same as in ( |75V Following the same procedure as that described at the beginning of 
IV Bl one finds that the perturbation of the power spectrum of v is given by the expression 

AV 

— = F cS = 25R[J (A) (109) 
ro 

with 

I(k) = 2i J dr,'^[v {0) (r,',k)] [(l-c 2 s )k 2 -a 2 m 2 cB ]v {0) , (110) 

where cj. and rri^g are explicitly given in (11041 i and d!05t , respectively. 

Using the Gaussian ansatz d89l for the turn, the deformation of the power spectrum is found to be (see Appendix for details) 

J 7 off (k) ~ Fo (k) + J~i (k) + O (™^ 4 ) , (111) 

where Tq is the contribution corresponding to the decoupling {rrih — > oo) limit. For large values of jj,, namely p, 3> x*, it is given 
by 

J r o(k) = - ^ ^ (x» sinx, + cosx») (sinx, — a;* cosx») . (112) 
The next order term F\, which is the main contribution due to the heavy mode, is given (again in the limit fx 3> a;*) by 

Fiik) = _ „ 5- (a;* sinx* + cosx*) (sinx* — re* cosx*) . (113) 

2v7r x% m% 

The behavior of the power spectrum is illustrated in Fig{T3]for different choices of parameters. For modes that exited the 
Hubble radius before the turn, i.e. x* -C 1, the power spectrum is almost constant. For large values of x*, corresponding to 
modes that exited the Hubble radius after the turn, the power spectrum oscillates around its unperturbed value, with the amplitude 
of these oscillations suppressed for large values of x*. The power spectrum reaches its maximum value around x* ps 1, i.e. for 
modes that crossed the Hubble radius during the turn. 

One also observes on the left panel that the power spectrum slightly changes with the value of rrih / H, in agreement with the 
above remark that the massive mode contribution is subdominant. An increase of mh/H corresponds to a slight decrease of the 
amplitude. By contrast, the power spectrum is much more sensitive on the parameter fl, as illustrated on the right panel. This 
means that the enhancement of the power spectrum is stronger for sharper turns, as expected. 

While we have focussed our discussion on sharp turns, i.e. ft ^> 1, it is worth mentioning that, in the opposite limit fl <C 1, the 
turn lasts for several e-folds and one recovers the constant-coupling case, discussed e.g. in In our notation, this corresponds 
to a constant 9 m w A9 \ij s/2n, leading to 

c 2 « l - 1 >f , m c 2 ff « - V > 1 - — . (114) 
Trm 2 2tt V m l J 

The resulting theory is analogous to a single field theory with slightly modified sound speed and spectrum tilt. 

Let us finally compare the results obtained from the effective theory with those derived from the two-field analysis. A crucial 
remark is that, in the decoupling limit rrih oo, the heavy mode contribution in (|99l becomes negligible and the result coincides 
with the dominant contribution (11 12l i of the effective theory. This confirms the expectation that when the turn is not too sharp, the 
effective theory is sufficient to describe the power spectrum. When the sharpness of the turn reaches a critical value, the effective 
theory breaks down and an additional contribution appears in the power spectrum. 
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FIG. 13: Power spectrum dl 1 U obtained from the effective theory. The plots of the left panel correspond to different values of 
rrih/H (with fixed parameters jl — 20 and A9 — 7r/10). The plots on the right panel correspond to different values of jl (with 
fixed parameters rrih/H — 100, A8 = 7r/10). 



VI. CONCLUSION 



In this work, we have studied the impact on inflation of heavy modes, which are expected to be ubiquitous in the context 
of high energy physics models of inflation that lead to a multi-field landscape. A crucial ingredient of our investigation is the 
use of the mass basis, which consists of the eigenvectors of the effective mass matrix for the linear perturbations, instead of the 
traditional kinematic basis. Within this basis, we have derived the multi-dimensional low energy effective theory governing the 
light degrees of freedom, starting from an initial system with an arbitrary number of light modes and of heavy modes. Our result 
thus generalizes the results of yflU], where a single light field was assumed. 

We have then explored the imprints of these massive modes when the inflationary trajectory undergoes a sudden turn. Assum- 
ing a simple Gaussian ansatz for the time evolution of the angular velocity, parametrized by the mass scale /i ~ 1/At tuin , we 
have been able to obtain an analytical description of the evolution of the velocity direction with respect to the lightest direction 
of the potential. The sharpness of the turn is governed by the ratio n/m^. When this ratio is of order one or higher, the turn is 
sharp and the angle between the velocity and the light direction undergoes oscillations, with frequency u ~ mn and damped by 
expansion (oc et~ 3 / 2 ), during the relaxation phase after the turn. We have confirmed this behavior by numerical integration for a 
concrete model. 

We have also investigated in detail how the turn affects the final adiabatic power spectrum. We have confirmed that, for a 
soft turn, the low energy effective description is sufficient to account for the power spectrum. However, for a sharp turn, the 
power spectrum contains several contributions, which we have been able to estimate analytically or numerically in our simplified 
description. The resulting power spectrum exhibits oscillations periodic in k and it would interesting to look for this type of 
signature in the CMB data, following previous investigations of oscillatory features (see e.g. ll25l - l28ll 5 ). 

Another aspect, which we have not considered here, is the impact of the tiny oscillatory component of the background due to 
the oscillating trajectory in case of a sharp turn. This effect, recently discussed in B0h32I1 . is complementary to our analysis. It 
would be interesting to estimate this effect with our formalism, taking into account explicitly the coupling between the light and 
heavy modes at the level of linear perturbations. 

Several extensions of the present work can be envisaged in the future. An immediate generalization would be to extend our 
results to non-standard kinetic terms, by introducing a non-trivial field space metric. It would be also interesting to refine our 
description of the turn by allowing additional light or heavy directions or by going beyond our simple ansatz for the angular 
velocity, in order to explore whether some features of the power spectrum are sensitive to the details of the turn. Finally, it would 
be worth going beyond the linear perturbations to investigate the specific imprints of heavy modes on non-Gaussianities. 
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5 Note that the reference t28il repeats the analysis introduced in f2^l , in the context of monodromy inflation, with WMAP7 data. 



21 



Appendix A: Analytical approximations for the power spectrum 

In this section, we give some details about the derivation of the analytical approximations given in the main body of the paper 
in Sec|Vl The goal is to solve analytically the integrals // and 1%, defined in ( f85l > and ( f86l >. by using the ansatz 

e m (t) = Ae-jLe- 1 ^ 2 * 2 . (ad 

V Z7T 



The correspondence between the cosmic time t and the conformal time r\ is given by 

a = -^- = e Ht =^> r) = rj*e- Ht (A2) 
where we fix = for simplicity. As a consequence, we can express B' m and 6" n as Gaussian functions of the cosmic time t: 

6' m (t) = a*A6-tj=e-^ 2t2+m , (A3) 

V Z7T 

C(t) = ^l»H (l - £t) e~i^ Ht . (A4) 
where a* = — (iJyy*) -1 is the scale factor at the turning point. 



1. Sub-horizons scales 



Deep inside the two horizons, the mode functions are approximated by 



u m « u m *- m = -j=- . (A5) 



Since the integrands in and Ih are exponentially suppressed when \t\ > /i , it is appropriate to use the approximation 

h 1 



e Ht « 1 — Ht + ^H 2 t 2 when i7 <C \i, which will be assumed here. We thus substitute in the integrals 



e ix,e- Ht e tx,(l-Ht+±H 2 t 2 ) 

u m « u h(0) « « 7 = , (A6) 



as well as 



* a: * ix,(l-Ht+^H 2 t 2 )+Ht I tt , zt2j\ e \n\ 

u h/ ^ a*— =e *v T 2 ;t {-H+H t) . (A7) 



4(0) 



/2fc 

The corresponding integrals /; and defined in fl85T l and ( |86t . can then be evaluated analytically. One finds 

(✓ \ _ l 4x 2 +,i 2 (l-<lx 2 ) x „ ( 1 - 4,i 2 +8/I 4 +4x J ) \ 

e^ + fi-ifi) e 4 (<^) +l , (A8) 

V M / / 

and 

40) = A#e^ (l-y&~ 2 ). (A9) 

2. Intermediate scales (1 < x* < i/) 

For the light mode, we start from the exact form of the mode function d75l l and use the same approximation for the exponential 
factor as in the previous subsection, which gives 

ikri f j \ e ix,(l~Ht+±H 2 t 2 ) 



22 



For the heavy mode, we use the approximation ((93), which yields 



_L fi -ln(-fcr ) ) v /q^ = _^^ e -lnx» e -(4+i,)Ht f ^ _ Vg^g 2 % < f -,ln 2 ) 



'2fc V2A; 

V2 J V2k 



One finds that the light integral can be written in the form 

-1/2 



h{k) 



i (A6>) 2 p 



with 



Ay/nx^ 

fi 2 -4: (fj 2 - 1) xj 

A{P 4 + x 2 ) ' 
a;* (1 - Ap 2 + 8/i 4 + 4x 2 
4(m 4 + x 2 ) 

And the heavy contribution is given by 

A8 



1 - ipj (e Pl+i9l x 2 + 2ie P2+ ' 1 " 2 x, - e P3+iq3 ) + e^ 1 (J 1 + x 2 ) 



Pi 



1>2 



P - (M " 2 ^ 



a* 4 + x'i 



P3 



9^ _ 4 ^ - 3 x* 



4(/2 4 + x 2 ) 



<12 



;i-2/i 2 + 2/2 4 +x 2 ) 
A 4 + 



93 



12/i 2 + 8/i 4 + 4x 2 



4(/i 4 + x 2 ) 



h(k) 



Ay/Vx 



3/2 



if (1— ln(2t>)) ^21/ In a;* 



\ -3/2 / \ -3/2 



with 



e Pi++i?i+ ( _ 3 _ 2» („ + a.) J + iei«+-Hfla+a;, U-^x* - 1 - 2i (f + a;„, , 



/_ (x*) = e? ,1 -+'"? 1 - 4-^ x* + 3 + 2i (v - x*) - *e 



x* ( 4— + 1 + 2i (y — x*) ] , 



where 



and 



Pi+ 

P2+ 
91+ = 
92+ = 

Pi- = 

P2- = 

91- = 

92- = 



9/i 2 — 4 (x* - 1/) [(A 2 - 3) - ^/i 2 

8(m 4 + x 2 ) 
p 2 — A (x* — v) [(p 2 -l) 



x* — i^/i 2 ] 



8(m 4 + x 2 ) 

-12^/i 2 + (-9 + 4j/ 2 + 12/i 2 - 8p 4 ) x* - 8«2 - 4xJ 
8(A 4 + * 2 ) : 
-4^/i 2 + (-1 + Av 2 + Ap 2 - 8p 4 ) x* - 8isxl - 4xJ 

8(A 4 + ^ 2 ) ' 

9p 2 - 4 (j/ + x*) {vp 2 + (p 2 - 3) x*) 
8(A 4 + x 2 ) ' 
p 2 - A (y + x*) (z^/i 2 + (/5 2 - l) x*) 

8(A 4 + ^) ' 
-12z//t 2 + (9 - Av 2 - 12/i 2 + 8p 4 ) x„ - 8i/x 2 + 4xJ 

8(A 4 + ^ 2 ) ' 
-4i^/i 2 + (1 - Av 2 - Ap 2 + 8p 4 ) x* - 8vxl + Axl 



8(p 4 + xl) 

In the limit p > x», (IA13b and (IA16b simplify into (J96j and (|97| i. 



(All) 
(A12) 

(A13) 

(A14) 
(A15) 

(A16) 

(A17) 
(A18) 

(A19) 
(A20) 
(A21) 
(A22) 

(A23) 
(A24) 
(A25) 
(A26) 
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3. Super-horizon scales (a* <S 1) 



In this case, d93l ) and thus ( IA1 II ) and ( IA121 > remain good approximations for the heavy mode. For the light mode, it is 
convenient to expand the expression for u/( ) with respect to the small parameter x« so that 

p -ik-q / • \ ix,e~ Ht / : \ 1 / A 1 



u m = —= 1 - — = 1 + — « -^= \ —e nt - -e-^ r < . (A27) 



After some manipulations, we get 

h = (Mf ^-j={i + ie^xl) , (A28) 

Ih = -M {2v+ ?\^ + < 1+ ^) x l'\^m. (A29) 
6\Jv 



4. Effective theory 



The initial condition for the mode function vim is the same as in dAlOb . namely 



V(0) (t, k) 



2k 



1 + — e 

a; 



HI 



(A30) 



Substituting in the integral I(k) defined in dl 10b , and using dA3BA4b with the expansions ( |106t and dl071 >, we find, at leading 
order, 



Io(k) 



(A9ff, 



tie 4 ! 1 + a^e 4 * 1 + 



and, at next-to-leading order, 

H 2 (A6) 4 fl 3 



\JP? - 13 



9i+4^3-2,l 2 + ii,)i 
, 4(iA 2 + x,) 



2 37 jj- C 



ia; 2 e 4i ^ 2 + 4x 



(A31) 



47T 3 / 2 TO 2 \/2/i 2 — iXirX 3 

l-4ix*+4a: 2 



2e s,i 2 -4,*„ ^ ^ie 2 ^ 
ff 2 (A0) 2 /2 



2e 4(2 l ^+x,) a;>t _ ix 2 



1 4^x„ / 1 \ 



(p? — ix*) 5 / 2 x 3 

: ^i e Th + T0f^ + T^ft (9 + 2/t 2 ) (/i 2 - ix*) 5/2 + 32e fj^ (ijl 



2 i \ 2 
1 



-16ie 4 ^- 4 ^*/ia; 4 (i/2 2 + a;*) 2 - ie 4 m 2 (^ 2 +*.) ^^2 _ ^^2 + ^ 2 (_15 + 2/t 2 + 16a; 2 ) 

l-8ix*+4x 2 

+e 4 a 2 - 4 -x, ^ (2i/i 6 + 52ia; 2 + 8i/2 2 a;* (13i + a;*) + /i 4 (-47i + 14a:* - 20«x 2 )) 

l-3ix t +x 2 

+4e ^ (_^ 6 - 18a; 2 - A/^x* (9i + a;*) + /i 4 (8 + 17^* + 10x 2 )) 

+ie ^-4^,*^ (_2/i 6 - 36.x 2 - 8/i 2 a;* (9t + a;*) + /i 4 (-9 + 54zx» + 20a; 2 )) }. 



(A32) 
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